Superconductivity near the vibrational mode instability in MgCNia. 
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To understand the role of electron-phonon interaction in superconducting MgCNia we have per- 
formed density functional based linear response calculations of its lattice dynamical properties. A 
large coupling constant A= 1.51 is predicted and contributing phonons are identified as displacements 
of Ni atoms towards octahedral interstitials of the perovskite lattice. Instabilities found for some 
vibrational modes emphasize the role of anharmonic effects in resolving experimental controversies. 
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The discovery of superconductivity in MgCNis has 
generated a new puzzle in the recent series of found 
superconductors [3- Despite its relatively low Tc'SK, 
the presence of Ni signals the possible importance of 
correlation effects which makes the physics of the pair- 
ing mechanism relevant to the famous high Tc cuprates 
and brings the discussion of unconventional non-electron- 
phonon mechanism. The experimental information char- 
acterizes MgCNia as moderate 0, or strong 0| cou- 
pling conventional superconductor through the analysis 
of specific-heat data, supports the s-wave pairing by 
CMR experiments [5|, and at the same time shows a 
zero-bias anomaly in tunnelling data Q . A clear need for 
detailed information about the phonon spectra emerges 
from these controversial data in order to clarify the role 
of electron-phonon interaction (EPI) and understand the 
mechanism of superconductivity. 

In this work, we perform theoretical studies of the 
strength of the electron-phonon coupling in MgCNia by 
using fully self-consistent density functional based linear 
response calculations of the lattice dynamical proper- 
ties as a function of phonon wavevector q. This method 
was proven to provide reliable estimates for the phonon 
spectra and electron-phonon interactions in a large vari- 
ety of systems 0,0- We extract a large coupling constant 
A = 1.51 and identify phonons contributing to it as dis- 
placements of Ni atoms towards octahedral interstitials 
of the perovskite-like cubic structure. We find some of 
the lattice vibrations to be unstable in linear order which 
emphasizes the role of strong anharmonic effects in the 
possible resolution of the recent experimental puzzles. 

The basic element of the perovskite MgCNis structure 
is given by a carbon atom placed at the center of the 
cube and octahedrally coordinated by 6 Ni atoms. Our 
electronic structure calculation using full potential lin- 
ear mufEn-tin orbital (LMTO) method 5| reveals Ni-d 
C-p hybridized valence bands ekj in accord with the pre- 
vious studies 0, 0, 0|. The Fermi surface consists of 
several sheets such as rounded cube sections centered at 
r, thin jungle gym area spanning from R [^H]^ to 
M [iiO] ^points, dimpled square shaped hole pockets 
centered around X point [^00]^ as well as little ovoids 
along r — _R.The tight-binding picture discussed before 
consists of carbon px Py Pz orbitals hybridized with 



d states of three Ni atoms numerated accordingly as Ni^; 
Niy and Ni^.For example (see Fig. 1) carbon state hy- 
bridizes with Ni'z dz2_i2 and also with Ni^; dxz Nfy dyz. 
Similar picture holds for carbon pa;,p„ orbitals which are 
not shown. As it has been noted this in particu- 
lar results in two antibonding states crossing the Fermi 
level, which in the nearest neighbor approximation have 
no dispersion along some directions in the Brillouin zone 
(BZ). 




FIG. 1: Basic element of the structure and set of tight-binding 
orbitals relevant to low-energy physics MgCNis. Arrows show 
displacements of Ni atoms corresponding to wavevector q = 

1-2 2"' a ■ 

The appearance of nearly flat areas of ekj gives rise 
to 2D van Hove singularity (vHS) placed 40 meV below 
the Fermi energy Ep which is responsible for a strong 
narrow peak near ep. This has generated speculation 
about closeness of MgCNisto ferromagnetic instability 
upon doping 0, [ijj, |ll|, |l^ . The narrowness of the vHS 
band is controlled by the second-nearest neighbor hop- 
ping integrals, which for the states shown in Fig.l corre- 
spond to the overlap between Ni^^ and Nij^ dxy orbitals. 
As we discuss in this work, nearly unstable phonon modes 
exists when each of the two Ni atoms moves toward oc- 
tahedral interstitial sites. (This is shown in Fig. 1 by 
arrows for a phonon wave vector q = [j^O]^). We shall 



2 



see that these distortions wipe out the narrow vHS peak 
and give rise to a large electron-phonon coupUng. 

To calculate lattice dynamics of MgCNia as a func- 
tion of wave vector q we utilize the linear response 
method^ I3, We use 2k LMTO basis set, generalized 
gradient approximation for exchange-correlation^^, ex- 
perimental lattice constant a=7.206 a.u., as well as ef- 
fective (40,40,40) grid in k-spacc (total 1771 irreducible 
k points) to generate the phonon dispersions LOqu and 
electron-phonon matrix elements ffk+qj'kjOn a (10,10,10) 
grid of the q vectors (total 56 irreducible q points) . 

Our calculated phonon spectrum along major high 
symmetry lines of the cubic Brillouin zone is given on 
Fig. 2. The frequencies are seen to be span up to 900 K, 
with some of the modes showing significant dispersion. In 
general, we distinguish three panels where the top three 
branches around 900 K are carbon based, the middle 
three branches around 600 K are Mg based and 9 lower 
branches are all Ni based. For the F point z-polarized 
modes, in particular, consist of (i) Ni^-C against Mg-Nij,- 
Nij, vibrations (186 K), (ii) pure Nia;-Nij, vibrations (262 
K), (iii) Mg against Ni^-Ni^ (561 K) (iv) C against Ni^ 
vibrations (857 K). 
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FIG. 2: Calculated phonon spectrum of MgCNia using density 
functional linear response method. Some curves are widened 
proportionally to the phonon linewidths. 

A striking feature of this phonon spectrum is the pres- 
ence of a low-frequency acoustic mode which is very soft 
and is even seen to be unstable along (^^0) direction 
in the BZ. This mode is essentially Ni based and corre- 
sponds to perpendicular movements of two Ni atoms to- 
wards octahedral interstitials of the perovskite structure. 
The latter is made of each of the four Ni atoms and two 
Mg atoms. For example, considering the xy plane (see 
Fig. 1) for the q point M such movements can be seen as 
a 2D breathing around this vacant interstitial. We find 



a similar situation for other wave vectors and in other 
directions of the BZ, where each pair of Ni atoms prefers 
such in-phase displacements perpendicular to each other. 
The softness and instability here can be understood as 
the octahedral interstitials are only places to escape for 
each Ni atom stressed between two surrounding carbons. 

The discussed displacements affect the overlap inte- 
grals between nearest Ni t2g orbitals which control the 
width of the vHS band. For the xy plane these are the 
hoppings between orbitals (see Fig. 1). It is therefore 
clear that these modes should have large EPI. As they 
bend the Ni-C bonds, this effect has already been pre- 
viously emphasized*^ by considering rotations of the Ni 
based octahedra at the R point. Here, we point out move- 
ments for wave vectors away from R point, where, for ex- 
ample, Nia;-C and Nij,-C bonds can be seen as changing 
their relative angle. 




FIG. 3: Calculated one-electron structure corresponding to 
the Ni based frozen phonon with q — [^^0]^. Top panel 
- undistorted bands, bottom panel - distroted bands corre- 
sponding to the Ni displacements by 0.2 A. 

To illustrate the crucial change in the electronic struc- 
ture due to such distortions. Fig. 3 shows two one- 
electron spectra in the vicinity of the Fermi level cor- 
responding to the M point frozen phonon involving Ni^. 
displacement along the y axis in phase with the Ni^ dis- 
placement along the x axis as illustrated on Fig. 1. The 
top panel of Fig. 3 corresponds to the undistorted en- 
ergy bands drawn in the original cubic BZ for easier 
comparison with the published data 0, llfil llll | . As the 
point [1^0]^ is now a reciprocal vector of the new dou- 
bled lattice, the bands are seen to be simply folded, and 
the narrow vHS band is readily recognized. When we 
introduce a distortion by 0.2 A, the only essential dif- 
ference is the width of the vHS band which now dis- 
perses as much as 1.5 eV. Given the smallness of the as- 
sumed displacement this emphasizes the large electron- 
phonon coupling. Its appearance cannot be understood 
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from a simple geometric overlap between the two t2g or- 
bitals {dxy states between Ni^j and Niy shown on Fig. 
1). We therefore look for an electronic enhancement 
due to nesting-like features of the Fermi surface. In- 
deed, such nesting can be found for the two dimpled 
square shaped hole pockets centered at X separated ex- 
actly by the wave vector [^^O]^. We have confirmed 
that feature by corresponding calculation of the integral 
^(^kj — eF)S{eu.+qj — cf) which provides the total 
phase space available for the electrons to scatter at given 
wave vector q with no energy change. 

We now turn our discussion to the detailed depen- 
dence of the electron-phonon coupling across the entire 
BZ. This at least can be done for all stable phonons. 
Fig. 2 shows the calculated phonon linewidths 7q;yby 
widening some representative dispersion curves ujq^ pro- 
portionally to 7q,y. Each phonon linewidth is propor- 

where the electron-phonon matrix element is found self- 
consistently from the linear response theory 0. In par- 
ticular, we see that some phonons have rather large 
linewidhts. This, for example holds, for all carbon based 
higher lying vibrational modes. The strength of the cou- 
pling, Aqi,, for each mode can be obtained by dividing 
7qy by 7rA^(ei?)a;qj,,where iV(ei?) is the density of states 
at the Fermi level equal to 5.3 st./[eV*cell] in our cal- 
culation. Due to large w^^, this unfortunately results in 
strongly suppressed coupling for all carbon modes which 
would favor high critical temperatures. The coupling, 
however, is relatively strong for the Ni based modes. For 
example, we can find A's of the order of 1-3 for the Ni 
based optical phonons around 250 K. Again, the anal- 
ysis of the polarization vectors shows that these vibra- 
tions involve Ni movements towards octahedral intersti- 
tials. For example, (see Fig. 1) the movement of Nia, 
and Nij, atoms along the z direction either in-phase or 
out-of phase result in larger overlap between Ni^ Nij, dyz 
orbitals and between Ni^ Ni^, d^z orbitals. Similar to 
what we find for the M point using the frozen phonon 
method (Fig. 3), this again enlarges the vHS bandwidth 
resulting in large EPI. A extremely large coupling (A^TO) 
occurs for the soft acoustic mode involving the interstitial 
breathing (M-point). Here we point out a triple effect: 
(i) breathing of four Nia, Nij, atoms into the interstitial, 
which results in larger dxy overlap, (ii) nesting enhance- 
ment which helps wiping out the vHS peak, and, (iii) 
the smallness of uj^^. As the linewidth of this particular 
phonon is so large, the concept of phonon itself has to 
be questioned, but due to a smallness of the phase space 
associated with this vibration, this has a little effect on 
integral characteristics such as A. 

Unfortunately, finding the integral value of A is another 
challenging problem due to the appearance of the imag- 
inary frequencies. Neglecting the unstable mode com- 
pletely results in the average coupling constant equal to 
0.95 mainly due to the discussed Ni vibrations around 



250 K (see Fig. 2).. It is however clear that the low- 
frequency mode has a large contribution to A and cannot 
be omitted. We have asked a question for which wave 
vectors q this mode is unstable? Since we know its dis- 
persion across the entire BZ, we can determine a sur- 
face in q space which separates the real and imaginary 
frequencies. Fig. 4 shows the result of such an analy- 
sis. We see the area around the F point which continues 
along the lines towards the R, M, and X points. Here we 
find the stability of the mode. The area around the M 
point, where we find enormously large EPI, is seen to be 
very small . This is also clear from the dispersion rela- 
tions shown on Fig. 2 for the (^^0) direction where the 
frequency becomes real just near the M point itself. 




FIG. 4: Surface in the Brillouin zone which separates stable 
and instable areas for the Ni based acoustic mode. 

This instability carefully avoids symmetry points 
which can be understood by keeping in mind the symme- 
try of the discussed distortions. Namely, anharmonicity 
is expected to be large for all in-plane movements in- 
volving two, three, or four Ni atoms towards octahedral 
interstitials. This is forbidden near the F point and for 
all wavevectors (^00) along FAl as well as for along 
TR directions if, for example, xy plane is considered. It is 
however allowed for (^^0) along FM and for {■^£,0) along 
XM. Moreover, for the wavevectors with ^ 0, the 
displacements away from xy plane are allowed and the 
instability is quickly suppressed. This, e.g., is the case of 
(iiC) along MR. 

The persistence of the instability which does not occur 
for any of the high-symmetry point needs a non-trivial 
frozen-phonon analysis. As our polarization vectors 
promt that the largest anharmonicity is expected for the 
TXM plane with ~ 0, we have performed three such 
calculations for the points a =(^0)^' ^ ^(HO)^' 
and M shown on Fig. 4. This corresponds to the dis- 
cussed displacements of two (point a), three (point 6), 
or four (point M) Ni atoms within xy-plane. The re- 
sults of these calculations reveal esentially anharmonic 
interatomic Ni potentials. A shallow double well with 
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a depth of the order of 20-40 K and the curvature at 
the equilibrium of the order of 40-50i K exists for the a 
and b points which becomes vanishing at the M point. 
We have found such a behavior by both total energy and 
force calculations for the supercells of 20 atoms (a and b 
points) and 10 atoms (M point). Clearly, such a small 
depth on the temperature scale indicates that the distor- 
tions are dynamical and zero point motions would pre- 
vent the appearance of the static long-range order. This 
picture is confirmed by our recent extended x-ray ab- 
sorbtion fine spectroscopy ex per iments (EXAFS) which 
reveal dynamical distortions |l5| not seen in the diffrac- 
tion experiments. Remarkably that the parameters of the 
double well (depth "'20K and the equilibrium curvature 
~50i K) extracted from EXAFS are in accord with our 
frozen phonon data. 

It is now clear that our instabile mode is not related to 
a statically distorted structure of MgCNis but should be 
resolved using an anharmonic theory of the EPI. Using 
a combination of our frozen phonon and linear response 
data we are able to map the behavior of the instable 
mode inside the entire Brillouin zone using a (4,4,4) grid 
of q-vectors. This gives us a possibility to determine a 
contribution to the Eliashberg spectral function using a 
generalized formula ll^ 

q kjk'j' n 

[/k, - h'rWl., - eko- +cc>„(q)]|Gg-,i,^.(q)P 

where the summation over the eigenstates n of the anhar- 
monic oscillator appears here with Ci;„(q) — e„(q) — eo(q) 
being the excitation frequencies around the ground state 
level for the anharmonic well at a given wave vector 
q. The generalized electron-phonon matrix element 
^k"]'kj('5) involves transitions from the ground to nth 
excited state and also includes changes in the effective 
one-electron potential to all orders with respect to the 
displacements. The most important is the effect of the 
frequency renormalization which appears in denomina- 
tor of Eq. As we find, our anharmonic mode has 
the first excitation frequency which, depending on the 
q-point, varies from 120 to 150 K. We also find that 
the summation over n is lastly convergent in accord with 
the previous work[l^: the oscillator strength for n=2 is 
always less than 2% of that for n=l. The linear electron- 
phonon scattering matrix elements are known to us from 
the linear response calculations, and numerically small 
corrections due to second-order couplings have been ex- 
tracted from our frozen phonon data using the band split- 
tings technique [l^ . 

Our resulting value of A for this anharmonic mode ap- 
pears to be 0.56. Adding the result for A = 0.95 from all 
other modes our total A is now 1.51. This is consistent 



with the values of A extracted from specific heat measure- 
ments and would cause MgCNia to be a strongly 
coupled electron-phonon superconductor. To estimate 
the Tc, we solve the Eliashberg equation with our total 
(harmonic-|- anharmonic) a'^F{uj). The cutoff frequency 
LUcut is taken to be 10 times the maximum phonon fre- 
quency. The Coulomb pseudopotential fi*{uJcut) is varied 
and the Tc within the range 7-20 K is obtained. To make 
connections with the values of fi* which are standardly 
discussed with the Allen-Dynes modified McMillan T^ 
expression ^| we need to use a rescaling formula 

= (t^cut)]'^^ + ln(wcut/a;iog) with wiog 120K. 
We found that that both the experimental Tc and the su- 
perconducting energy gap A(0) can be reproduced with 
the value of /i* = 0.33. While the effect of spin fluctu- 
ations on superconductinvity and mass enchancement is 
better to discuss in terms of its own coupling constant 
Aspm, at the absence of detailed calculation of the lat- 
ter, we can only conclude that our enchanced /i* can in 
principle appear due to localized nature of Ni d orbitals. 
It is, for example, known that spin fluctuatlions are im- 
portant in V and fi*~0.3 is needed'2^. In fact, at the ab- 
sence of spin fluctuations, fi* is usually 0.1-0.15 and Tc in 
MgCNiawould be larger. We conclude that spin fluctua- 
tions partially suppress superconductivity, the result ex- 
pected from the conventional theory. Note that the same 
conclusion has been reached based on the recent analysis 
of the specific heat data 0|. We finally calculated the 
phonon contribution to the specific heat and fitted it to 
the form C{T) — j3T^ at low temperatures. Our value 
of /3 = 0.35 mJ/[mol*K'*] is close to the values 0.40-0.42 
mJ/[mol*K'*] determined experimentally j3l ij. IT^ . 

In conclusion, by performing linear response calcula- 
tions of the electron-phonon interaction in MgCNia we 
reported the value of A = 1.51 consistent with the strong- 
coupling limit of electron-phonon mechanism of super- 
conductivity. The unusually large anharmonic correction 
to A for the lattice near instability is emphasized. 
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